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ABSTRACT 

Delayed detonations of Chandrasekhar-mass white dwarfs (WDs) have been very suc- 
cessful in explaining the spectra, light curves, and the width-luminosity relation of 
spectroscopically normal Type la supernovae (SNe la) . The ignition of the thermonu- 
clear deflagration flame at the end of the convective carbon "simmering" phase in the 
core of the WD is still not well understood and much about the ignition kernel dis- 
tribution remains unknown. Furthermore, the central density at the time of ignition 
depends on the still uncertain screened carbon fusion reaction rates, the accretion his- 
tory and cooling time of the progenitor, and the composition. We present the results 
of twelve high-resolution three-dimensional delayed detonation SN la explosion simu- 
lations that employ a new criterion to trigger the deflagration to detonation transition 
(DDT). The simulations fall into into three ignition categories: relatively bright SNe 
with 5 ignition kernels and a weak deflagration phase (three different central densi- 
ties), relatively dim SNe with 1600 ignition kernels and a strong deflagration phase 
(three different central densities) and intermediate SNe with 200 ignition kernels (six 
different central densities). All simulations trigger our DDT criterion and the resulting 
delayed detonations unbind the star. We find a trend of increasing iron group element 
(IGE) production with increasing central density for all three categories. The total 
^^Ni yield, however, remains more or less constant, even though increased electron 
captures at high density result in a decreasing ^^Ni mass fraction of the IGE mate- 
rial. We attribute this to an approximate balance of ^^Ni producing and destroying 
effects. The deflagrations that were ignited at higher density initially have a faster 
growth rate of subgrid-scale turbulence. Hence, the effective flame speed increases 
faster, which triggers the DDT criterion earlier, at a time when the central density 
of the expanded star is higher. This leads to an overall increase of IGE production, 
which off-sets the percental reduction of ^^Ni due to neutronization. 
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1 INTRODUCTION 

SNe la have come to fame as the Universe's most lumi- 
nous standardizable candles - crucial i ngredients to the 
study of dark energ y and cosmology (e.g. iRiesjfet al.lQQ^ : 
ISchmidt et al.|[l998l ). A limiting factor on the precision of 
using SNe la as distance indicator s is the inherent scatter i n 
their normalized light curves (e.g. IWood-Vasev et al.ll2007l ). 
A better understanding of the intrinsic vari ation of super- 
nova brightnesses and spectra is needed (e.g. lAlbrecht et al.l 
l2006l : iMiknaitij 120071 '). Simulations of SN la explosions 
are already being used to aid in improving the preci- 
sion of cosmological distanc e measurements bas ed on su- 
pernovae in the future (e.g. iBlondin et al.|[201ll ). In addi- 



tion, SNe la also play a critical rol e in galaxy gas kine- 
mati cs (e.g. Sca nnapic co et al. [2008^ . positron production 
(e.g. Chan fc Li ngcnfc ltei) llQQa ). and chemical evolution 
(e.g. iMatteucci fc Greggid [l98^ ) . Detailed modeling of the 
explosions is therefore useful for understanding the ori- 
gin of the Galactic 511 keV line, the origin and evolu- 
tion of heavy elements, and kinetic supernova feedback, 
and measuring the Hubble parameter as a function of red- 
shift. The standard model of SNe la relies on the nu- 
clear fusion of the initial composition (predominantly ^^C 
and ^"^O) of a massive white dwarf (W D) star to more 
tight ly bound nuclei to power the explosion (|Hovle fc Fowleij 
1 19601) . The exact nature of the progenitor systems and de- 
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tails of the dynamics of the nuclear burning processes how- 
ever are not known. Among the leading scenarios are the 
Chandrasekhar-mass models, in which a WD accretes mat- 
ter from a companion star and grows in mass to near the 
Chan drasekhar limit u ntil pycnonuclear carbon fusion reac- 
tions (jCameronl Il959l ) start taking place. Once carbon fu- 
sion reactions produce more energy than is carried away 
by neutrino losses, the core becomes convective and when 
the nuclear burning time of a fluid element becomes shorter 
than the eddy turnover time a deflagration flame may be 
born fe.g. lWooslevlll990l ). Numerical simulations of the con- 
vective stage leading up to the ignit io n of th e deflagration 
were p erformed by iHoflich fc SteinI (I2OO2I). iKuhlen et al.l 
ll2006l).|Piro fc Chand (|2008h . IPiro fc BildstenI (|2008l ). and 
IZingale et al.l (|2009l ). The central density of the WD de- 
creases significantly during the simmering phase between the 
onse t of carbon burning and the ignition of the deflagration 
(e.g. iLesaffre et al ] |2006l : IPiro fc Bildstenll2008l l. Those cal- 
culations, however, are not taking electron captures and the 
URCA process correctly into account, and some uncertainty 
in the evolution remains. The rate of the s creened ^^C-^^C 
fusion reaction is still quite uncertain (e.g. Itoh et al.ll20oj : 
IJiang et al.l l2007l : ICasgues et al] |2005| . |2007| ). The "central 
density at the time of ignition, however, depends only mildly 
on the exact value of this r eaction rate (|Cooper et al.|[2009l : 
llapichino fc Lesaffrd I2OI0I ). More important is the initial 
mass and the accretion and cooling history of the WD, which 
determines the thermodynamic state of the interior. This re- 
sults in a range of possible central dens ities at ignition, from 
less than 2 x lO'' to over 5 x 10^ g cm"^ (jLesaffre et al.ll2006h . 
Metallicity has a considerable impact on the supernova 
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shown to depend rather w eakly on metallicity and the CO 
ratio iLesaffre et aLlliood ). If the ini tial deflagration flame 
can transition into a detonation (e.g. Khokhlov et al.lll997l : 
lRopkell2007l : IWooslevllioOTi : IWooslev et al.ll2009l l. then good 
agree ment of the models with o bservations can be obtained 
(e.g. 'Ropkc fc Ni emeveJ l2007l : iBravo fc Garcia-Senj I2OO8I : 
iKascn ct al. 2009). A successful explosion model has to re- 
produce the observed range of peak absolute magnitudes 
(i.e. ^®Ni masses) and the width-luminosity relation and 
scatter thereabout. Furthermore, the observed correlation 
between the brightness of an event and the delay time or 
age of the host stella r population has to be explained (e.g. 
[Gallagher et al.ll2008l \ Recently, a connection between age 
of the host stellar population and SN la brightness was pro- 
posed vi a the effect of longer cooling times on the ignition 
density (jKrueger et al.ll2010l ). Varying the central density 
for 150 two dimensional delayed detonation supernova sim- 
ul ations within th e statis tical ignition framework presented 
in iTownslev et all (j2009l ). the authors found that the '^^Ni 
yield decreased with increasing central density, while the to- 
tal iron group element (IGE) yield remains roughly constant. 
This is attributed to increased production of stable isotopes 
(such as e.g. ^''Fe or ^*Ni) due to increased neutronization 
via electron captures at the higher densities. There are, how- 
ever, at least three competing effects that influence the '''^Ni 
mass produced in a delayed detonation SN. 

(i) Electron capture rates on protons and iron-group iso- 
topes under electron degenerate condisitons are strongly 



incre asing with density (e.g. iLanganke fc Martfnez-Pinedd 
I2OOII) . Consequently, a distribution of nuclei in nuclear sta- 
tistical equilibrium at high density neutro nizes at a much 
faster rate than one at lower density (e.g. ISeitenzahl et al] 
l2009t ). which acts to lower the ^^Ni mass. 

(ii) Near Chandrasekhar-mass WDs in hydrostatic equi- 
librium with a higher central density are more compact, i.e. 
signiflcantly smaller and slightly more massive and tighly 
bound. This may translate into a more compact WD at the 
time of the flrst DDT, which could lead to an overall larger 
part of the WD being burned to IGEs, which acts to raise 
the ^"^Ni mass. 

(iii) Deflagrations e volve differently at higher gravita - 
tional acceleration g (|Khokhlovl Il995l : I Zhang et alll2007l ). 
From linear stability analysis, the Rayleigh- Taylor temporal 
growth rate scales with y/g. The different flame evolution 
and turbulence generation could have an effect on the DDT 
(e.g. the transition density), which, depending on the dif- 
ferent degree of "pre-expansion" , could either lower or raise 
the ^'^Ni mass. 

The effect of variations in the central density of the WD 
on pure deflagr ations has been exp lored in three-dimensional 
models before (|Ropke et al.ll2006l ). Here, we present the re- 
sults of twelve high-resolution three-dimensional delayed det- 
onation SN la simulations (that employ a new DDT crite- 
rion, see Section 12. 3p for three different ignition conflgura- 
tions and a range of central densities. We flnd that, for the 
same spatial ignition spark distributions, the ^^Ni yield re- 
mains more or less constant a function of central density at 
ignition. The deflagrations that were ignited at higher den- 
sity produce subgrid-scale turbulence at a higher rate, which 
triggers the DDT criterion earlier when the central density 
of the star is higher. This leads to an overall increase of 
IGE production as well as enhanced electron captures. Even 
though the mass in ^^Ni comprises a smaller fraction of the 
mass that has burned to IGEs, the overall ^®Ni yield remains 
roughly constant since more total mass in IGEs is produced 
in the detonation. Only the cases where much of the IGEs 
are produced in the deflagration phase show a trend of de- 
creasing *®Ni with central density. In Section[2]we introduce 
our setup and briefly review the computational methods, in 
Section |3] and |4] we present and discuss the results, and in 
Section [5] we conclude. 



2 METHODS AND SIMULATIONS 

The large computational demands of the high-resolution 
three- dimensional simulations we perform prevented a 
statistical framework approach similar to th e one pre- 
sented in iTownslev et"aLl (|2009l ) . [Kriieger et al.1 l|2010h and 
I Jackson et al. I (|2010l ). Under the constraints of limited com- 
putational resources, we chose six different densities for a 
setup with an intermediate number of ignition points (200 
kernels), and three densities each for the setups with the 
least (5 kernels) and the most (1600 kernels) ignition points 
respectively. The central densities are such that they cover 
the distribution of ignition densities e xpected from differen t 
cooling ages and accretion histories ILesaffre et al.l (|2006h . 
The ignition spark configurations are selected in a way that 
SNe with a range of brightnesses with ^®Ni masses between 
~ 0.45 and 1.1 M0 are obtained. 
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2.1 Initial models 



2.2 Computational method 



All simulations presented here are full star simulations per- 
formed in 3D. The initial stellar models are cold, isothermal 
(r = 5 X 10^ K) WDs in hydrostatic equilibrium with cen- 
tral density Pc ranging from 1.0 to 5.5 x 10® g cm^'\ The 
composition is assumed to be 47.5 per cent ^^C, 50 per cent 
^^O, and 2.5 per cent ^^Ne (to account for solar metallicity 
of the zero-age main-sequence progenitor) by mass homoge- 
neously throughout the star, resulting in an electron fraction 
Ye = 0.49886. 



A strong deflagration phase leads to more energy re- 
lease and hence expansion of the star. The ensuing deto- 
nation then produces less ^^Ni, leading to a dimmer event. 
In the multi-spot ignition scenarios, the strongest deflagra- 
tions are obtained by placing an optimal number of igni- 
tion sparks approximately s ymmetrically about the center 



tion sparks approxima.tely s ymmetrically about tne center 
jGarcia-Senz fc Bravj l2005l : iLivne et aiTl2005l : iRopke et all 
l2007l ). While too few ignition sparks lead to an overall weak 
deflagration, too many of them lead to vigorous burning in 
the initial stage and thus an early expansion of the WD 
that suppresses bu rning in later stages of the deflagration 
jRopke et al.ll2006l ). 



Asymmetric ignition spark distributions lead to a 
weaker d eflagratio n phase and henc e a brighter SN la 
ijRopke et al. .20071 : iKasen et al.ll2009l ). For the ignition of 
the deflagration we use setups generated from a Monte- 
Carlo based algorithm. The primary input parameters are 
the number of the ignition kernels and the distribution 
type. The details of the ig nition process rema i n unk nown. 
IWooslev et al] (|2004l ') and IWunsch fc WooslevI (|2004l ') con- 
clude from analytical models that multi-spot ignition within 
the inner ~150 km or so is a possible scenario. The total 
number and spatial distribution of the ignition spots, how- 
ever, was not conclusively constrained by their models. We 
investigate three different sets of explosion models corre- 
sponding to different ignition scenarios. We choose config- 
urations of 5, 200, and 1600 kernels which are spherically 
arranged around the center of the WD following a Gaus- 
sian distribution in radius. The placement of kernels with a 
distance greater than 2.5 times than a given variance a is 
suppressed. For the setups with 5, 200, and 1600 kernels, we 
set a — 0.6, 0.75, and 1.8 x 10^ cm, respectively. The ra- 
dius of the spherical ignition kernels is set to i?k = 10® cm. 
Finally, we impose a length scale Dk, which the distances 
between the centers of the ignition kernels have to exceed. 
Dk is set to 10® , 3 x 10^ , and 5 x 10* cm for the setups with 5, 
200, and 1600 kernels, respectively. Note that for i?k > i'k, 
the sparks may partially overlap which is the case for the se- 
tups with 200 and 1600 kernels. Within a given model suite 
(i.e. 5, 200, or 1600 kernels), the locations of the ignition 
sparks are only once randomly determined in the beginning; 
the resulting spatial realization of the ignition configuration 
is then kept fixed and identical in all the simulations with 
different central density. With these choices of ignition spark 
distributions we cover a large range of 5®Ni masses, between 
~ 0. 45 and 1.1 M^, which is consistent with norm al SNe la 
(e.g. IContardo et"al]|2000l : IStritzinger et al.ll2006al lbl). 



The reactive Euler equations are solved using a finite volume 
scheme based on the PROMETHEUS code bv lFrvxell et all 
(jl989f l. which is an imple mentation of the "p i ecewis e 
parabolic method" (PPM) of IColella fc WoodwardI (|l984h . 
The grid resolution is 512 x 512 x 512 cells for all sim- 
ulation s. We use the expan d ing hy brid grid implementa- 
tion of iRopke fc Hillebrandtl (|2005al lbl) . with a uniform in- 
ner grid that contains the deflagration level set and a non- 
uniform outer grid that covers the remainder of the com- 
putational domain. Our simulation code is based on a large 
eddy simulation (LES) approach, which resolves the largest 
turbulent structures and models the turbulence on unre- 
solved scal es using a turbulence subgrid-scale model (for 
details see ISchmidt et~all l2006allB . The code uses a co- 
movine grid (|Ropkell2005l : iRopke et al.l[2006h with an outer 
coarse grid following the WD's expansion and an inner finer 
grid tracking the flame front. The flame itself is treated as 
a discontinuity separating fuel and ash; its propagation is 
tracked with the level se t technique dOsher fc Sethian|[l988l : 
ISmilianovski et al. 1 1 19971 : iReinecke et al.lll999i r In this thin 
flame approximation, the energy liberated in the nuclear 
burning is released immediately behind the level set rep- 
resenting the flame surface. Since nuclear matter burned 
in a deflagration undergoes different burning than matter 
processed in a detonation, separate leve l -set representations 
are u sed (|Golombek fc Niemeveij l2005l : iRopke fc Niemevej 
|2003). Using a full nuclear reaction network in every com- 
putional cell to calculate the source terms for the hydro- 
dynamics is currently still computationally too expensive 
for three-dimensional simulations. We solve this problem by 
tabulating the energy release as a function of f uel density. 
For the detonation we use the new tables from iFink et al.l 
(|2010l ). A table for the energy release of the deflagration 
level set was calculated in a similar way. 



2.3 DDT criterion 

The transition from a subsonic deflagration to a super- 
sonic d etonation based on th e Zel'dovich gradient mech- 
anism (IZel'dovich et al.l Il970l) was introduced to SN la 
theory b y iBlinnikoy fc Khokhlwl (119861) and further ana- 
lyzed by [K hokhlovl (119913 Ibl): iKhokhlov et all ([l993) and 
iNiemever fc WoosIct] (|l997h . The main result of their stud- 
ies was that such a transition is only possible in the tur- 
bulent deflagration stage, where large velocity fluctuations 
v' lead to a mixing of cold fuel and hot ash up to a cer- 
tain length scale. These "hot sp o ts" ar e supposed to be the 
seeds of a DDT. iLisewski etlll (|2000l ') pointed out that v' 
must exceed 10* cm s~^. Indeed, velocity fluctuations on 
this scale have already been found in three-dimensional de- 
flagration simulations ([Ropkc 2007). Wooslcy ( 2007 ) argued 
that for DDTs there are specific restrictions on the burning 
properties deep in the distributed burning regime, which is 
the regime where strong turbulent flame interactions are ex- 
pected. As DDTs cannot be resolved in full-star simulations, 
we employ a subgrid-scale model to calculate the probability 
of these transitions. 

The details of this subgrid-scale DDT model, which is 
guided by the late st studies of the micr oscopic mechanism 
of DDTs in SNe la lWooslev et all (|2009l ). are described in a 
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separate paper (Ciaraldi-Schoolmann & Ropke, in prepara- 
tion). It accounts for the intensity of the turbulent velocity 
fluctuations as well as the fuel density pfuoi and fuel frac- 
tion Xfuei in the grid cells crossed by the flame front. If 
the probability P{v' > v'^-^^^) to find velocity fluctuations 
larger then 1)^^;^ in a specific area ylflamo at the flame front 
exceeds a certain threshold Acrit, detonations are ignited 
in the grid cells which contain the largest velocity fluc- 
tuations. Aflamo is defined as the part of the flame where 
Pfuci = [0.6 - 0.8] ■ lO'' g cm"^ and Xfud = [0.3 - 0.7]. To 
properly estimate this area, we take the fractal dimension of 
the flame front into account, which is ~ 2.36 ([Sreenivasanl 
Il99ll : iKersteinI [l988l : I WooslevI [ioOTi ') . The number of igni- 
tions is given by the ratio Aflame to Acrit- The criterion must 
hold at least for half of an eddy turn Tcddyiya ~ ^/^(-^)) 
where L is the turbulent integral scale a,nd V(L ) the veloc- 
ity at this scale. Following conside ratio ns oflWooslcy ( 2007[ ) , 
ICiaraldi-Schoolmann et all (|2009l ) and lRopkel (|2007h we as- 



sume L = 10*^ cm and V{L) « [10^-10**] cm s"\ We choose 
a constant va lue of Teddy, /o = 0.00 5 s i n our analysis . We 
further follow iLisewski et all (|2000l ') and iRopkel (120071 ) and 
deflne D^rit ~ 10* cm s~^ and Acth = lO^'^cm^ as our thresh- 
olds for the DDT-criterion. While the details of the imple- 
mentation are beyond the scope of this publication and will 
be presented elsewhere, we point out that this modeling ap- 
proach is significantly different from simply fixing a certain 
DDT threshold density. Our criterion in addition requires 
strong local turbulent velocities. It is thus more restrictive 
and substantially reduces the number of DDTs. 



3 RESULTS 

The chosen distinct setups lead to different evolutions of 
the deflagration flame. In turn, the different evolutions of 
the flame front has great impact of the properties of the 
following delayed detonations (see Fig. [T]). 

On the one hand, the energy released during the de- 
flagration stage up to the time tooT when the first DDT 
occurs, is smaller for simulations starting with fewer igni- 
tion kernels (see Table [1]). The ensuing weaker expansion 
manifests itself in a higher central density a,t t = turn- This 
in turn translates into a larger fraction of the total mass of 
the WD burned into IGE during the detonation phase (see 
Table[2} than for the models with more ignition kernels. We 
thereby confirm that the strength of the defiagration is a pri- 
mary parameter for the ^®Ni production (|Ropke fc Niemeved 
I2OO7I ). and hence brightness, of a supernova that explodes 
in the delayed detonation scenario (see Fig. [2]). A strong 
deflagration results in fainter events, primarily due to the 
large expansion of the star prior to the DDT, but also due 
to the copious neutronization in the defiagration phase (see 
Fig.©. 

On the other hand, for an identical spatial distribu- 
tion of ignition kernels, the simulations with higher initial 
central density exhibit stronger turbulence production and 
the subgrid-scale energy grows at a faster rate initially (see 
Fig. [4|. This can be understood by the difference in gravi- 
tational acceleration g. For the same spatial distribution of 
ignition spots, the distribution in the mass coordinate will 
be centered further out at larger mass and thus larger g in 
the higher central density case. For Rayleigh- Taylor dom- 




3.0 X 10' cm 



f =teDT = 0.679 s 




10* cm 



Figure 1. Shown arc snapshots of the deflagration level set 
(blue) and the detonation level sot (red) for the lowest (left col- 
umn, pc = 1.0 X 10^ g cm"'') and the highest (right column, 
Pc = 5.5 X 10^ g cm"'') central density, a) Both simulations ini- 
tially have a spatially identical arrangement of the 200 ignition 
kernels, b/c) After t = 0.3s, the deflagration has burned signifi- 
cantly more for the high density case, d/e) The spots where the 
DDT criterion is first triggered are circled; the detonation triggers 
at an earlier time for the high density case, f/g) The detonation 
level set is propagating through unburned fuel away from the 
DDT spots. Multiple detonations can be launched as long as the 
DDT criterion is fulfilled. In spite of the differences in time evo- 
lution and morphology, both models produce the same amount of 
56Ni. 



inated deflagration flames the tur bulent burning velocity 
scales with (e.g. lKhokhlov|[T995l l. Consequently, at equal 
time after ignition (e.g. 0.3s), a high central density model 
will have burned signiflcantly more mass than a low central 
density model (compare Figs.[l]3 and[l]:). As a further conse- 
quence, the higher the initial central density, the sooner our 
DDT criterion is fulfilled (see Table [1] but note the outlier 
with 5 ignition kernels for pc ~ 3.5 x 10^ g cm~^). Impor- 
tantly, there is a trend that models with high initial central 
density also have a higher central density at t = tom, which 
is a proxy for the amount of fuel at densities high enough 
that it will be burned to IGE in the detonation. We find 
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Table 1. DDT attributes for all models. Tabulated are the nuclear energy released -Enuc, the central density pc, and the (average) density 
of the first DDT spot(s) f5i(t = fDDx) at the time (ddt when the first DDT(s) occured. 



Model 


Pc(i = 0) 


iDDT 


Enuc{t = iDDT) 


Pc(i = iDDx) 


Pl(i = iDDT) 




[10^ g cm-3] 


[s] 


[10^1 ergs ] 


[10* g cm-3] 


[10'' g cm-3] 


0005 


1.5 


1.253 


0.334 


4.911 


0.773 




3.5 


0.890 


0.246 


11.610 


0.758 




5.5 


0.911 


0.488 


6.183 


0.716 


0200 


1.0 


1.135 


0.587 


1.450 


0.700 




1.5 


0.993 


0.667 


1.706 


0.746 




2.9 


0.802 


0.803 


2.270 


0.756 




3.5 


0.756 


0.838 


2.533 


0.749 




4.0 


0.755 


0.936 


2.193 


0.759 




5.5 


0.679 


0.930 


2.711 


0.761 


1600 


1.5 


1.077 


0.813 


0.718 


0.705 




3.5 


0.848 


0.990 


0.827 


0.779 




5.5 


0.757 


1.087 


0.875 


0.755 



Table 2. Nucleosynthetic yields for all models. Tabulated are the total WD mass Mtot, and the final masses of l^c, l^o, intermediate 
mass elements, iron group elements, and ^^Ni {Mi2Q,Mi6Q,MiMEyMlGE) a-nd Afse^i)- Furthermore tabulated are the masses of iron 
group elements and ^^Ni at the time iDDT when the first DDT(s) occured (Mjq^ and Af^g"^.), as well as their respective relative fractions 

'"lilkand*^-' 



of the final masses, ( 



A/lGE 



Model 


Pc(i = 0) 

[10^ g cm-3] 


A/tot 
[Mq] 


Ml2c 

[Mq] 


A'fisQ 
[Mq] 


A/lME 

[M©] 


MiGE 

[Me] 


A'/se Ni 
[Mq] 


[Mq] 


[Mq] 


A-flGE 


. rdcf 


0005 


1.5 
3.5 
5.5 


1.378 
1.406 
1.416 


0.004 
0.002 
0.002 


0.033 
0.017 
0.015 


0.205 
0.100 
0.098 


1.136 
1.287 
1.301 


1.033 
1.139 
1.095 


0.221 
0.168 
0.288 


0.182 
0.125 
0.172 


0.19 
0.13 
0.22 


0.176 
0.110 
0.157 



0200 



1.0 
1.5 
2.9 
3.5 
4.0 
5.5 



1.361 
1.378 
1.400 
1.406 
1.409 
1.416 



0.017 
0.012 
0.008 
0.006 
0.007 
0.007 



0.104 
0.073 
0.063 
0.053 
0.057 
0.053 



0.447 
0.390 
0.346 
0.307 
0.314 
0.280 



0.793 
0.904 
0.984 
1.040 
1.031 
1.076 



0.698 
0.768 
0.752 
0.770 
0.735 
0.698 



0.385 
0.447 
0.547 
0.573 
0.622 
0.626 



0.321 
0.348 
0.359 
0.354 
0.375 
0.325 



0.49 
0.49 
0.56 
0.55 
0.60 
0.58 



0.460 
0.453 
0.477 
0.460 
0.510 
0.466 



1600 



1.5 
3.5 

5.5 



1.378 0.014 
1.406 0.015 
1.416 0.015 



0.104 
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Figure 2. Shown is the mass of '''"Ni produced in the different ex- 
plosions as a function of the total nuclear energy liberated during 
the deflagration phase up to i = iDDT- 



that for all ignition kernel distributions, the total yield of 
IGE material within an ignition distribution suite increases 
with Pc (see Fig.[5]and Table[2|. The total yield of ^"^Ni ap- 
pears flat with Pc for the model suites with 5 and 200 ignition 
kernels; only the model suite with the strongest deflagration 
phase (1600 ignition kernels) has a trend of decreasing ^®Ni 
with Pc (see Fig. [g] and Table 12]). 



4 DISCUSSION 

We attribute the almost constant *'^Ni mass to a coinciden- 
tal balance of the competing effects presented in Section [TJ 
- the larger electron capture rates at higher central den- 
sity are offset by a greater total mass in IGE due to the 
inherent compactness of the WD and faster evolution of the 
flame towards DDT. The high density simulations with 5 
and 200 ignition sparks exhibit more subgrid-scale energy 
generation at early times and therefore higher flame speeds 
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Figure 3. Shown is the relative mass fraction of ^^Ni to IGE pro- 
duced in the different explosions as a function of the total nuclear 
energy liberated during the deflagration phase up to t = tDDT- It 
is evident that the strength of the deflagration (as measured by 
En_uc{t = ^ddt)) is a very good proxy for the mass ratio of ^®Ni 
to IGE. 




Figure 4. Shown is the subgrid-scale turbulent energy as a func- 
tion of time for the model sequence with 200 ignition kernels. Note 
that the rate of turbulent energy production is initially larger for 
the high density cases. 



Figure 5. Shown is the mass of IGE as a function of central 
density pc of the WD at the time the deflagration was ignited. 
A trend with increasing IGE mass with central density is evident 
for all three ignition configurations. 
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Figure 6. Shown is the mass of ^^Ni as a function of pc{t = 0). 
Evidently, for these simulations the ignition configuration is the 
primary parameter that determines the ^®Ni mass (and hence 
peak brightness) of the events, whereas the central density is 
merely a secondary parameter. 



initially. The still highly turbulent deflagration flame reaches 
the outer layers of the WD with low fuel density faster, and, 
consequently, the DDT criterion triggers earlier when the 
central density of the star is still higher (see Table [T]). As a 
result, these simulations produce more IGE and even though 
the ^®Ni fraction of the IGEs is lower due to increased neu- 
tronization, (see Fig. (H)), the total amount of ^^Ni remains 
roughly constant. The difference between ^^Ni and IGE is 
mainly made up of stable iron group nuclides such as ^*^Ni, 
^''Fe, but also other radioactive nuclides such as ^^Ni or ^^Co 
contribute. 

The model suite with 1600 ignition kernels, which has 
such a high density of ignition sparks that the whole central 
region is filled with burning products of the deflagration, be- 
haves differently. Due to the numerous ignition sites, a large 
part of central mass of the star is already burned in the de- 
flagration before the first DDT occurs. The IGE produced 
in the deflagration, where most of the electron captures oc- 



cur, are a large fraction of the total IGE produced (see Ta- 
ble [2]). The electron captures and resulting shift of ''^Ni to- 
ward more neutron rich stable Fe-group isotopes occur most 
copiously behind the slowly moving deflagration flame front, 
and consequently the increased production of IGE material 
does not reflect in larger ^''Ni masses for cases where the 
deflagration contributes most of the IGE mass. The strong 
deflagration and vigorous expansion leads to such low cen- 
tral densities at t = tDDT that the ensuing detonation cannot 
produce sufficient ^^Ni to counter this trend. Events produc- 
ing such large amounts of strongly neutronized IGE matter 
cannot, however, make up most SN la events, due to the 
unusual isotopic composition (e.g. IWooslevlll997l '). 

We can only speculate why our simulations predict 
increasing IGE and roughly constan t ^^Ni product i on for 
higher initial central density, whereas iKrueger et al] (|201G| ) 
find the opposite - decreasing ^^Ni and constant IGE. 
One possible reason for the different trends between the 
two sets of simulations is the way the DDT is handled. 



SN la diversity from central ignition density variations 7 



For simulations of delayed detonations in SNe occurring 
via a DDT, the detonation is Keneral l y put in "by hand" 
re.g.lArnett fc Livne"l994': 'Livnc"l9 99l:lGamezo et al.ll2005l : 
iBravo fc Garcia- Scnz. 2008 : Krucgcr et al.ll2010l ) . usually by 
choosing a critical densit y where a deflagratio n transitions 
to a detonation. Recently, I Jackson etaD (|2010l ) investigated 
the effect the particular choice of such a transition density 
has. They found a quadratic dependence of the IGE yield 
on the log of the transition density. Instead of imposing a 
fixed transition density, we utilize a dynamic DDT crite- 
rion (see Section 12. 3p . which takes the effects of different 
deflagration evolutions on the detonation initiation into ac- 
count. We point out that typical densities where our DDT 
criterion trigg e rs (se e Table [T]) are lower than 10^ g cm~"^. 
Ijackson et al.l ()2010h have shown that the variance of the 
^^Ni yield for a statistical set of simulations is relatively 
large for such a low choice of transition density (see figure 3 
from their work), in agreement with our observed large range 
of 5®Ni masses obtained. Numerous other obvious differences 
between the simulation sets exist, including the nature of the 
propagation and the nuclear energy release of the burning 
fronts (level sets vs. reaction progress variables) , the dimen- 
sionality of the simulations (3D vs 2D), or the structure of 

the computational mesh (AMR vs. ex panding grid). 

In this context, note the work of lMeakin et al. I (|2009D . 
who present a suite of supernova explosion models with dif- 
ferent off-sets for the initial defiagration bubble. Although 
their single bubble off center ignition scenario does not ex- 
plore central density at the time of the ignition of the defia- 
gration as a parameter, their result that a strong defiagra- 
tion phase need not necessarily result in less ^^Ni produced 
is the same. They also find that the total amount of IGE 
decreases for models that had a more vigorous defiagration 
phase (leading to more expansion) before the detonation is 
triggered, but the ^^Ni yields remains approximately con- 
stant (see figure 12 of their work). 



5 CONCLUSIONS 

We have performed twelve three-dimensional hydrodynam- 
ical simulations for delayed detonation SNe la for a range 
of central densities and ignition conditions. We find a trend 
of increasing IGE production with central density within 
each set of ignition conditions. This is because the high cen- 
tral density WDs are more compact and the fiame evolves 
faster; the DDT occurs sooner when more unburned ma- 
terial is still above the density threshold (« 10^ g cm^'*) 
where a detonation will still produce IGE. In spite of the 
larger IGE mass, the more vigorous neutronization occur- 
ring in the high density models during the defiagration phase 
yields ^^Ni massest hat are more or less constant with pc for 
the brighter SNe. Only dim SNe, which have a strong de- 
fiagration phase and expansion prior to the DDT, exhibit a 
trend of decreasing ^^Ni mass with increasing density, since 
the increased neutronization in the defiagration phase can- 
not be compensated for by the relatively weak detonation 
phase. This trend, however, is of secondary importance when 
compared to the effects of varying the ignition kernel dis- 
tribution. For a given ignition kernel spatial distribution, 
the central density therefore infiuences the brightness of the 
supernova event only as a secondary parameter. From the 



works of iTownslev erall (|2009l ') and iBravo et al.l (|2010l ), it 
appears that the same holds for composition, i.e. metallicity 
and C/O ratio. Indeed, based on an analysis of high-quality 
V and B-band li ght curves of SNe la from the Carnegie Su- 
pernova Froiect. lHHflich et al.1 (|201G| ') propose that the com- 
position and central density are two independent secondary 
parameters for SN la light curves. In light of the impor- 
tance of the ignition configuration of the defiagration for the 
brightness of the SN, it is most crucial to establish how the 
central density at ignition (cooling time) and metallicity af- 
fect the statistical properties (notably number and location) 
of the ignition sparks themselves, and not their respective 
direct effects on the outcome of an explosion once a random 
ignition spark distribution was chosen. One should there- 
fore aim to quantify which effect composition, cooling and 
accretion history have on the ignition process, for example 
by mapping them into the exponenti ation parameter Cc of 
the sto chastic ignition prescription of ISchmidt fc Niemeved 
(|2006j) . This would require a better understanding of the 
physics leading up to ignition, including the nature of the 
convection and effects of electron captures and the convec- 
tive URCA process. 
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